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NASA-GSFC IONOSPHERIC CORRECTIONS TO 
SATELLITE TRACKING DATA 


ABSTRACT 

The purpose of this report is to present an overview of 
the development, verification, and recent implementation of 
the NASA-GSFC ionospheric model for satellite tracking data 
corrections. This highly successful model has been incor- 
porated into the Goddard Trajectory Determination System 
which is providing continuous trajectory computation support 
for the lunar orbiting Radio Astronomy Explorer-B launched 
on 10 June 1973. 
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NASA-GSFC IONOSPHERIC CORRECTIONS TO 
SATELLITE TRACKING DATA 


SUMMARY 

The Earth's ionosphere introduces systematic perturbations to all radio 
tracking data. At present NASA radio tracking is conducted in the approximate 
frequency range of 140 MHz to 6 GHz. At the lower frequencies daytime iono- 
spheric biases can easily reach several kilometers in range, several meters 
per second in range rate, and up to two or three milli radians in angle. At the 
higher frequencies these values decrease as the inverse square of frequency. 
However in all cases this ionospheric bias is greater than the basic tracking 
system resolutions involved (e.g. ref. 1). The lower tracking frequencies are 
preferred from the standpoint of simplified antenna pointing and hence ease of 
signal acquisition. The higher frequencies are preferred from the standpoint 
of high telemetry data rates and extremely precise orbit computation such as 
associated with satellite geodesy. 

This report describes the general development and salient features of the 
highly successful ionospheric correction model incorporated into the recently 
completed Goddard Trajectory Analysis Program (GTDS) which eventually will 
provide operational support for all GSFC missions. The Radio Astronomy 
Explorer-B (RAE-B) launch in June 1973 was the first of the series of NASA 
missions to be supported by GTDS. RAE-B is being tracked by the VHF (136 - 
148 MHz) Goddard Range and Range Rate System and ionospheric corrections to 
all tracking measurements are being automatically applied on a continuous basis. 
These corrections are a function of predicted solar activity, tracking station 
latitude and longitude, universal time and station to spacecraft geometry. 
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NASA-GSFC IONOSPHERIC CORRECTIONS TO 
SATELLITE TRACKING DATA 


1.0 EVOLUTION OF IONOSPHERIC RESEARCH 

Over the past fifty years thousands of reports and technical papers have 
been compiled regarding radio wave propagation phenomena as influenced by 
the Earth's ionosphere. The "non-exhaustive" bibliography (cited as Ref. 2) 
lists 1404 references for the period 1923 to 1960 alone. In addition basic iono- 
spheric theoretical considerations have been documented in hundreds of excellent 
books such as represented by references 3 through 5. 

Since the early days of radio (Circa 1925) ionospheric research has had 
practical application in terms of skywave radio propagation (Ref. 6). The effects 
of the ionosphere on radio astronomy, which had its beginning around 1930, have 
also been of long standing interest. The advent of radio tracking of artificial 
satellites in 1957 added still another motivation for a more complete understand- 
ing and prediction of ionospheric phenomena. As an aside it is also of interest 
to note that ionization in the upper atmosphere is directly linked to solar activity 
and solar activity has been monitored in terms of Zurich sunspot number since 
1755! Over the intervening 218 years a total of twenty consecutive eleven year 
sunspot cycles have been observed (Ref. 7). 

In early 1970 it was decided by NASA-GSFC that this wealth of theoretical 
and experimental background should be applied toward the development of a 
worldwide ionospheric correction model consistent with anticipated tracking 
accuracy requirements of the 1973 and beyond time frame. The philosophy was 
to start with the promising closed form mathematical profile description first 
suggested by Kazantsev (Ref. 8) and later modified by Bent in 1967 (Ref. 9). 

This profile was then to be related to empirically determined electron density 
decay factors derived from analysis of many thousands of satellite and ground 
based soundings, and satellite electron density measurements. The data of 
interest was already available for a seven year time span and stored in a number 
of data center archives. The ready availability of suitable large scale high speed 
computers in 1970 gave assurance that such a goal involving large scale data 
analysis was indeed feasible. By mid 1971 the basic NASA-GSFC ionospheric 
model had been developed and an extensive testing and evaluation took place. 
Since integrated electron content is the scaling parameter for ionospheric 
corrections to all satellite tracking data, every means available was used to 
compare this model with measured integrated content. The latter measurements 
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included: Faraday rotation from ATS-1 and ATS-3, group-phase velocity 
measurements using NASA Explorer series spacecraft, two-frequency ranging 
to ATS-3, intercomparisons with the U. S. Army Map Service SECOR System, 
and ray tracing through composite profiles derived from bottom and topside 
sounding measurements. In other tests corrections were applied to numerous 
passes of previously uncorrected tracking data and postflight orbit computations 
were performed with the general result of a significant reduction in residuals. 

RMS observed minus calculated range, range rate or angle data were often re- 
duced by factors of 4 or more. In addition many trajectory overlap tests were 
run which further indicated significant improvements in orbit computation accu- 
racy. These tests included comparisons with existing but much cruder ionospheric 
correction models and here accuracy improvements by factors of two or more were 
observed (Ref. 10, 11). The final step was merging of the total NASA-GSFC 
ionospheric model with the Goddard Trajectory Determination System (GTDS) in 
time for successful operational support of the Radio Astronomy Explorer-B 
launched in June of 1973. The time scale from the initial concept of this project 
to final implementation was approximately 3 years. As implied in the initial 
comments, this work necessarily reflects an extension and application of the 
efforts of many researchers. However this particular NASA contribution is 
directly associated with those persons cited in Figure 1 and in particular with 
Dr. R. B. Bent and Ms. S. K. Llewellyn who were the principal persons involved 
in the analysis and formulations which led to the present NASA-GSFC ionospheric 
model. This work represents a part of the ongoing applied atmospheric physics 
effort (Figure 2) in the Geodynamics Program Division as directed toward such 
areas as high precision satellite radio tracking and satellite altimetry in support 
of such programs as the NASA-GSFC Earth and Ocean Physics Applications Program. 


2.0 DEVELOPMENT OF THE NASA-GSFC IONOSPHERIC MODEL 

2.1 DEVELOPMENT OVERVIEW 

The types and quantity of measured data involved in the development of the 
NASA-GSFC ionospheric model are indicated in Figure 3. The starting point 
for the f q F 2 and M(3000) F 2 worldwide time dependent spherical harmonic 
mapping was already available as a result of many years of research at the 
former National Bureau of Standards, Boulder, Colorado , now the National Oceanic 
and Atmospheric Administration - Boulder (Ref. 3 and Ref. 12 through 19) . 

As will be explained, f G F 2 is the vertical sou ndin g frequency related to 
the ionosphere profile maximum electron density. The M(3000) F 2 scalar factor 
(Figure 3) is the ratio of the MUF(3000) F 2 to f o F 2 frequencies. The MUF 
(3000) F 2 is the highest frequency for a 3000 km single hop Earth point-to-point 
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NASA-GSFC IONOSPHERE MODEL DEVELOPMENT 



Figure 3 
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radio communication via the F 2 layer and is associated in this model with the 
height of the maximum electron density. 

Finally the world map of critical frequency f Q F 2 which is directly related 
to maximum electron density is updated with predicted (or measured if available) 
solar flux activity on a daily basis from, for example, the National Space Science 
Data Center at NASA-GSFC. The basic coefficients of the spherical harmonic 
expansion of the world map are stored in the computer memory and called up as 
required. The values computed from these basic stored coefficients are then 
updated as a function of measured or predicted solar activity. Figure 4 indicates 
this overall basic implementation. 


2.2 MODEL PARAMETERS 

The mathematically closed form ionospheric profile is shown in Figure 5. 

It consists of segments of parabolic, parabolic-squared and exponential functions. 
The functions are matched in magnitude and slope at each boundary. The maxi- 
mum electron density, N m , as will be shown, is related to the f o F 2 frequency. 

This density is on the order of 10 n to 10 12 electrons per meter 3 depending on 
time of year, time of day, solar activity, geographic latitude and longitude and 
so on. The corresponding vertical total integrated electron content varies from 
10 16 to 10 18 electrons per meter 2 . Also the height of the maximum density, 
h m , is a function of the scalar factor M(3000) F 2 and will typically range from 
200 to 400 km. The parameters affecting profile shape (e.g. Kj, K 2 , K 3 , Y m , Y { ) 
were empirically derived from the large data base indicated in Figure 3 and are 
stored in the computer in the form of look up values which are functions of lati- 
tude, longitude, f 0 F 2 and time of year. 

At this point a few brief comments regarding the nomenclature, origin and 
meaning of f Q F 2 and M(3000) F 2 is appropriate. It is well known by the com- 
munications engineer working with radio frequencies below 30 MHz that radio 
signals are reflected by apparently discrete ionospheric layers commonly termed 
E, sporadic-E, Fj, F 2 and so on. The daytime D layer is primarily an absorption 
layer (Ref. 20) that for example prevents broadcast signals (500 to 1600 kHz) from 
"Skipping out" to distant points during daylight hours. The profile in Figure 5 ob- 
viously does not model these discrete layers; however, the maximum density 
corresponds to the most important reflecting layer, namely the F 2 region. 

Since satellite tracking biases result from integrated effects through the entire 
ionosphere the lower ionosphere fine structure can be generally ignored. The 
origin of the nomenclature D, E, or F layer is worth mentioning as a matter of 
historical interest. In the early experimental work of Sir Edward Appleton his 
first ionospheric reflections were obtained from an apparent layer which he 
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NASA-GSFC IONOSPHERE MODEL IMPLEMENTATION 



Figure 4 
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simply called the E layer corresponding to the notation of the reflected E vector 
of the electromagnetic radiation. Later in the winter of 1925 he obtained reflec- 
tions from two other layers below and above the E layer height and termed them 
respectively D and F. He reasoned there might be still additional layers both 
above and below the originally discovered E layer and hence retained this nom- 
enclature rather then renaming the layers, for example A, B, C and so on (Ref. 21). 
Later the daytime F layer was found to split into two layers which are now commonly 
termed Fj and F 2 . The F 2 layer corresponds to the layer of maximum iono- 
spheric electron density. The f x F 2 frequency, or often termed MUF (ZERO) 

F 2 , is the highest vertical sounding frequency which will be returned to Earth. 

The relationship between f o F 2 and f x F 2 is described in section 2.3. The nom- 
enclature is again traced to radio communications applications, where zero 
means vertical or zero skip distance and MUF is "maximum usuable frequency" 
for reflection. Similarily the MUF (3000)F 2 is the maximum usuable frequency 
for a communications single hop of 3000 km from the F 2 layer. The latter is 
directly related to the height, h m , of the peak of the F 2 layer. 

The computer stored modeled parameters thus fall into two categories. 

These are the profile shaping factors determined empirically and the factors 
which describe the worldwide maximum density and height of maximum density, 
primarily f D F 2 and M(3000)F 2 as updated with either predicted or measured 
solar activity. Where: M(3000)F 2 =MUF(3000)F 2 + f D F 2 . 

Solar activity is obtained from two data bases - namely, Zurich sunspot number 
and solar flux as measured at the 10.7cm wavelength. The sunspot history, as 
indicated previously goes back to the year 1755, and is derived from an observer's 
numerical count of sunspots. The two data types have been shown to reflect the 
same information and one empirical relationship (Ref. 22) between the two is 
given by: 


0 = (63.75 + 0.728 R + a 00089 R 2 ) 10“ 22 


watts/meter 2 /Hz 


( 1 ) 


Where: 0 = 10.7 k flux 

R = average Zurich sunspot number 
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2.3 MAXIMUM ELECTRON DENSITY AND HEIGHT OF MAXIMUM DENSITY 


The familiar Appleton - Hartree equation relating index of refraction in a 
charged particle plasma in the presence of a magnetic field is traced back to at 
least 1927 (Ref. 23). This equation has been applied to ionospheric phenomena 
in a variety of approximate forms and the approximation validity is largely de- 
pendent on the particular portion of the radio frequency spectrum of interest. For 
example, at frequencies of 140 MHz neglecting the effect of the Earth's magnetic 
field introduces a maximum error of 2% and this error decreases rapidly for 
higher frequencies. At typical vertical sounding frequencies (f 0 F 2 ) ranging 
from 5 MHz to 15 MHz the absorption terms can be ignored for refraction cal- 
culations but the Earth's magnetic field must be included. Reflection for vertical 
sounding occurs whenever the index of refraction tends to zero and this can be 
calculated for a given density from the simplified Appleton - Hartree equation 
given below: 


8- 


= 1 - 


1-i 

■ 2 


\ Y 2 /(l-X)± |y}/(1-X) 2 + Y 2 


1/2 


Ionospheric index of refraction 
squared neglecting absorption 

( 2 ) 


c = free-space velocity of light 
e = charge on an electron 
f = oj/2tt = frequency of a radio wave 
H = geomagnetic field intensity 
m = mass of an electron 
N = number density of free electrons 
e 0 = electric permittivity of free space 1 
6 = angle between wave normal and geomagnetic field 
A. = 2ttc/u> = free-space wavelength of a radio wave 
ix. = real part of the refractive index 
Hq = magnetic permittivity of free space 
to = angular frequency of a radio wave 
< 4 , = yu. 0 He/m = angular gyro frequency of an electron 
co L = aJjj cos 8 

a> N = yNe 2 /e Q m = angular plasma frequency 

coj. = co H sin 8 
X = co^/co 2 
Y = co H /co 

Y l = <w l/ £0 
Y T = CjO^/cO 

References: 23, 24 
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For purposes of calculating maximum electron density at the F 2 layer 
the conditions leading to /x 2 = 0 in equation (2) must be examined at the f 0 F 2 
frequency. The presence of the static geomagnetic field, makes the ionosphere 
anisotropic or doubly-refracting, for radio waves. A radio wave of general 
polarization, on entering the ionosphere, will be split up into two waves each 
having a characteristic polarization. The relative amplitudes of the waves will 
depend on the original polarization. Each wave will then travel with its own 
characteristic velocity and attenuation. The two waves will, in general, have 
elliptic polarizations in opposite senses and are termed the "ordinary" and 
"extra-ordinary" waves. Each wave is reflected at a level where its refractive 
index tends to zero in the absence of collisions. Hence each wave will have its 
own 'critical' or 'penetration' frequency which is the maximum frequency which 
can be reflected by the layer. 

The index of refraction for the two waves is given by taking the sign before 
the square root in equation (2) as positive and negative for the ordinary and 
extra-ordinary respectively. The ordinary critical frequency f 0 F 2 and the 
extra-ordinary critical frequency f x F 2 are related by the formula 

W = < f x F 2> ! — VW (3) 


where f H is the gyrofrequency of electrons. 

f H = — | B | = 2.84 x 10 10 I B | H 
H 2-nm 


(B in webers/meter 2 ) 


when: 


f 0^ 2 > > 

f x F 2»f H 
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we can rewrite equation 3 as: 


V 2 ~ *0*2 " 0 


(4) 


For calculating maximum electron density f F 2 is usually used since it is 
most easily measured. For calculating skip communications via the F 2 layer 
fF 2 is usually mapped, using equation 4 to obtain £F 2 from measured fF 
(Ref. 17). Figure 6 is a world map of the gyro frequency which can be used to 
convert between the f o F 2 and f x F 2 frequencies. 

The expression for maximum electron density is obtained from equation 2 
at the ordinary critical frequency (i.e. f Q F 2 ) as: 


N 

max 


(f„ F 2) 2 

80.6 


electrons/meter 3 


(5) 


(where f 0 F 2 is in Hz) 


The height of the maximum density is calculated on the basis of the equations 
of Appleton and Beynon (Ref. 25). The required M(3000)F 2 or M factor (also 
see section 2.2) predictions are obtained from a generalized data set issued by 
the National Oceanic and Atmospheric Administration. The equation used for the 
height of maximum electron density in the NASA-GSFC model is given by: 


where 


h = 1346.92 - 526. 40M + 59. 825M 2 km 

m 

MUF( 3000) F, 

M 

f 0 F 2 


( 6 ) 
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2.4 SPHERICAL HARMONIC MAPPING OF f/ 2 AND M(3000)F 2 

As indicated in Figure 4 a basic unchanging set of harmonic coefficients are 
stored in the computer for both f Q F 2 and M(3000)F . These coefficients reflect 
both a time and geographical dependence. The f 0 F 2 set is scaled by a quadratic 
function of solar flux. The M(3000)F 2 is scaled by a linear function of sunspot 
number. The choice of update parameters is a matter of historical development 
and the entire update information can in fact be tied to solar flux measurements 
if deemed desirable. In addition, as a result of the analysis of the 400,000 world- 
wide bottomside ionograms, a latitude dependent bias in f 0 F 2 was discovered 
relative to previously presented data by the National Oceanic and Atmospheric 
Administration (NOAA- Boulder, Colorado). Thus this modelling effort which for 
f F 2 and M(3000)F depended heavily on the earlier work at Boulder (Ref. 12 
to 19), has provided a marked refinement for those geographical regions signif- 
icantly remote from the U. S. 

For both the f 0 F 2 and M(3000)F 2 the form of the basic harmonic expansion 
is: 


K 

A(0A,T) = ^] D k( T > G k (tf.M 

k = 0 


( 7 ) 


where: 4 > (Geographic Latitude) 

\ (Longitude) 

T (Time) = local mean hour angle -180°< T< 180° 

The time dependent Fourier Series modulating the geographical term 
models the diurnal and seasonal variations and has the form: 

H 

D k (T) = A( 0 k > + (A^cos.T + B ( j k) sin j T) (8) 

j = i 


where 


H = 6 for f c F 2 
H = 4 for M(3000)F 2 


The latitude and longitude spherical harmonic terms are expanded as as- 
sociated Legendre functions much in the same manner as the geopotential function 
used in all orbit computation. 
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The number of terms used in this harmonic representation is given in 
Table 1. The number of terms selected takes into account the maximum number 
warranted in light of the data base quality to which these functions have been 
fitted. 


Table 1 

Harmonic Coefficients 


Number of Terms 


Time 

Series 

Geographic 

Terms 

Total 

f F„ 
0 2 

13 

76 

988 

M(3000)F 2 

9 

49 

441 


Figure 7 presents the type of map resulting from such a representation. 
This particular example is an f x F 2 monthly mean prediction for August 1963 
at UT 0600 hours. With the present NASA capability, daily updates for f Q F 2 
and M(3000)F 2 can be made and this was actually done for the critical first few 
days following the June 1973 launch of the Radio Astronomy Explorer-B. 


2.5 INDEX OF REFRACTION FOR FREQUENCIES ABOVE 100 MHz 

The satellite tracking data acquired at frequencies above 100 MHz will be 
biased while traversing the ionosphere because the index of refraction is other 
than that of free space (i.e. /x = 1.0). Equation 2 for frequencies above 100 MHz 
can then be approximated by: 


M = 



(9) 
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where: N e 


e 

m 


CO 


Electron density (electrons/meter 3 ) 
Electron charge = 1.602 x 10“ 19 coulombs 
Electron mass = 9.11 xlO" 31 kilograms 
Free space dielectric constant 
8.855 x 10" 12 farads/meter 
27rf (radians/sec) 


The tracking data distortion experienced relative to free space, is related 
to the integral of the electron density along the ray path. In all cases the sys- 
tematic distortion will be scaled by the integrated content (electrons/meter 2 ) 
along the ray path, however the exact manner in which the radio measurements 
are affected will depend on the details of the tracking system involved. 


2.6 INTEGRATED ELECTRON CONTENT PREDICTION 

The profile shown previously in Figure 5 is easily integrated in closed form 
and thus permits calculation of predicted total vertical integrated content over 
any geographical region. In satellite data corrections however the total integrated 
content of interest is that along the path between tracking station and satellite. 

It has been found that a good estimate of the total path delay can be established 
by using the vertical integrated content divided by the sine of the local elevation 
angle, E*, at the point the ray path intercepts the maximum electron density. 

The geometry involved is straightforward and with reference to Figure 8, is as 
follows : 


E* 


arccos 


o 

cos E 

a + h 

m 


( 10 ) 


where: a 


h 


m 


E 


E* 


Earth radius 

height of maximum electron density 
elevation angle at tracking station 
local elevation angle at 


Since Doppler bias results from the time variation of the index of refraction 
another parameter used is the elevation angle rate at the height of the maximum 
which is given by: 
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Figure 8 


where: E = elevation angle rate as observed at tracking station 

E* = elevation angle rate at h m 

The NASA model thus computes the intercept of the tracking station to 
satellite ray path with the maximum density layer in the ionosphere for any 
station location on Earth and provides the appropriate integrated electron 
content which is in turn used in all subsequent tracking data correction equations. 


3.0 IONOSPHERE MODEL VERIFICATION AND IMPLEMENTATION 

3.1 TRACKING DATA CORRECTION EQUATIONS 

As indicated earlier, the exact form of the ionosphere tracking data correc- 
tion equations is a function of the specific tracking systems employed. However, 
certain simplifications are possible in order to gain insight into the problem. 

The results derived in the following are in fact applicable to NASA range and 
range-rate tracking as a consequence of the manner in which phase and frequency 
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relationships are maintained throughout the satellite tracking link. The following 
will indicate mathematically how range and range rate biases are introduced and 
the results will be shown to be consistent with the often presented formulation 
for range bias: 


AR = 


40.3 


f 2 


N e dR meters 


(12) 


f = frequency in Hz 

N g = electron density along ray path electrons /meter 3 

integral taken along ray path between tracking station and satellite 
In the model described: 


S 


N dR = 


sin E ! 


,* e 1 ec t rons /me t e r 2 


(13) 


(where I y is verticle electron content 
calculated at ray path intercept with 
region of maximum electron density) 

All NASA range and range rate systems employ narrowband phase modulated 
CW systems having the form: 


E = A cos (a> c t + /3 sin V) 


(14) 


o) c = 2 7if c = carrier angular frequency 
- 2-n-f m = modulation angular frequency 
/3 = modulation index 

The modulation index, /3, is typically maintained between 0.5 and 1.5 radians. 
The NASA tracking systems employ transmission from ground station to a space- 
craft transponder which in turn retransmits a radio signal to the ground station 
(e.g. Ref. 26). Thus the ionosphere is traversed twice and two frequencies are 
involved since the transponded frequency is offset from the spacecraft received 
frequency. The two frequency effect can in general be considered by using: 


19 




1 


1 


( 15 ) 


f 2 




down 


However the precise formulation for the effective frequency, f, must be made 
in light of the specific tracking system considered (Ref. 27). 

For illustrative purposes the frequency translation can be ignored and the 
radio signal demodulated as if the ionosphere were traversed twice by a coher- 
ently transponded carrier signal at frequency, f . 

The narrowband phase modulated signal can be considered to consist of 3 
principal frequencies — namely, a carrier and two symmetrically situated side- 
bands (Ref. 28). The form of such a signal with amplitudes normalized for sim- 
plicity is: 


E = cos co t + cos (co + co ) t 

o' v o 


(16) 


- cos (co - CO 

v o m / 

As seen by equation 8 the index of refraction is frequency dependent and the 
effect on the range measurement and coherent Doppler phase can be calculated 
with reference to Figure 9. The functional representation includes a phase de- 
tector for range measurement determination and a Doppler extractor for range 
rate determination. In principal the demodulation requires the multiplication of 
signals by means of a non-linear device and undesired higher order terms are 
filtered out. Anyone familiar with the actual mechanization of a system such as 
the Goddard Range and Range Rate System will realize the gross oversimplifica- 
tion depicted by Figure 9. However the mathematics of the more complex system 
can be shown to reduce to an equivalent coherent spacecraft turnaround with a 
unity transponder frequency translation ratio (Ref. 26). Referring back to Figure 
9, when the indicated multiplications are performed and the desired signals 
extracted one obtains : 
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and the two-way integrated electron content is modeled as: 
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It should be pointed out that equations 17 and 18 are again approximations 
(for example Ref. 29) which neglect higher order terms. However the error 
introduced is negligible at the frequencies considered being less than 0.3% in 
the worst case. 

The interesting fact to note with regard to the ionospheric effect in equations 
17 and 18 is that whereas the range tone measurement experiences a phase delay 
the carrier phase measurement undergoes an effective phase advance. The range 
tone delay is used as a measure of the range to the spacecraft. The carrier phase 
shift is differenced over a short time span (typically 0.5 seconds) to provide an 
average Doppler frequency which is directly proportional to range rate. It is 
apparent that the effect of the ionosphere on the two types of measurements is 
of opposite sign and hence should in principle allow one to solve for integrated 
electron content by comparing range differences with integrated Doppler and/or 
differentiated range measurements directly with range rate measurements. 

Such comparisons have in fact been made and used as one means of evaluating 
the NASA-GSFC Ionospheric model. The difficult task has been selection of 
appropriate data smoothing, editing, and compression techniques. This is because 
the inherently coarser range tone (typically 20 kHz at VHF) phase measurement is 
to be compared with the high resolution carrier phase data which can, with 5 second 
averaging, measure range rate to a resolution of 0.5 cm/ sec at 140 MHz and to better 
than 0.1 cm/sec at 2 GHz (Ref. 1). 
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While angle data from antenna pointing is also available from range and 
range -rate sites it is seldom used in NASA orbit computation simply because 
the data quality is not consistent with the high accuracy range and range-rate 
measurements. Often range -rate only or range only orbit solutions are performed 
depending on satellite geometry and distance of spacecraft from the tracking site. 

The NASA 136 MHz Minitrack interferometer measures direction cosines 
and for this system ionospheric corrections are applied in terms of angle of 
arrival. Unfortunately no simple formulation exists to relate this ionospheric 
correction. It will suffice to say that the ionospheric ray bending (in contrast to 
tropospheric effects) is predominant at 136 MHz and this is discussed in reference 
30. The formulation adopted for the Goddard Trajectory Determination System 
to calculate angular deviations is presented in Ref. 31. Ionospheric effects 
in satellite-to-satellite tracking are discussed in reference 32. 


3.2 GROUP-PHASE PHENOMENA 

The results presented as equations 17 and 18 can also be arrived at con- 
sidering group and phase velocity effects through a dispersive media. Using 
such an approach leads to the conclusion that these equations also apply, again 
to a certain approximation, to pulse systems and amplitude modulated systems. 

It is beyond the scope of this report to present a detailed discussion of such 
concepts as group velocity and phase velocity. However further discussion 
is presented in Ref. 3, 33, and 34. 

The important consideration herein is experimental. That is, the obser- 
vability of such phenomena as predicted by equations 17 and 18. 

Figure 10 is obtained from a comparison of range and range-rate data and depicts 
extraction of the ionosphere induced range rate bias while using the NASA-GSFC 
136-148 MHz range and range-rate tracking system. As the elevation angle scan 
rate tends to zero (at the relative time of 150 seconds) the range rate bias also tends 
to zero. This is the result predicted by the time derivative of equation 18 if all 
changes in total electron content along the propagation path can be attributed to the 
elevation angle scan. This is a reasonable assumption for near Earth spacecraft, 
however for deep space corrections, where the maximum scan rate is essentially 
Earth rotation of 7.3 x 10“ 5 radians/sec, the time rate of electron content change 
due to diurnal effects must also be included. 

Figure 11 shows the group-phase comparison for a 500 second data stretch 
where data was recorded at a 4 per second rate. The data compression utilizing 
polynomial smoothing compressed the data over 10 separate sub-arcs and the 
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OBSERVED RANGE RATE IONOSPHERE BIAS (CM/SEC) 


IONOSPHERIC BIAS AT VHF TRACKING FREQUENCIES 

SANTIAGO 22 JANUARY 1972 VHF GRARR 
14 h 57 m 40 s (START) 
m 15 h 01 m 47 s (STOP) 



Figure 10 
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Figure 11 



resultant curve is termed "observed". Experience based on both simulated and 
observed data indicates that near Earth polynomial smoothing of data must be 
applied with caution to prevent errors introduced by the smoothing process itself 
(e.g. Ref. 35). The VHF range-rate data has been noted to be consistently of 
higher quality than range data in the sense that compressed range rate values 
time tagged at the same instant are generally independent of the degree of fit 
once polynomials over a cubic are employed and data stretches are limited to 50 
seconds or less. Range data on the otherhand is to a certain extent sensitive to 
degree of fit and the modulation of "observed" values in Figure 9 is attributed to 
low frequency cyclic noise on the range data. This range noise appears to be a 
function of spacecraft spin and range tracking loop performance under high 
dynamics (i.e. perigee passes of near Earth satellites). However it must be 
mentioned that the measurements being performed are far beyond the initial 
design intent for tracking system performance. In fact if range data could be 
processed to achieve the same accuracy for range rate estimates as the carrier 
determined range-rate there would be no further need for range-rate measure- 
ments. 

Figure 11 also indicates the importance of considering horizontal electron 
density gradients in modeling of the ionosphere for the purpose of correcting 
tracking data bias. In this particular example the predicted gradient shown in 
Figure 12 accounts for approximately 20% of the predicted range rate bias shown 
in Figure 11. The effect of horizontal gradients is further demonstrated by 
Figure 13. In this case a typical simulated satellite overhead pass is considered 
and the range bias at 140 MHz compared for polar and equatorial passes. The 
difference between the two results again indicates the importance of modeling 
horizontal gradients. 


3.3 TWO-FREQUENCY RANGING 

Another well established means for experimentally determining integrated 
electron content is to perform satellite range measurements at two widely 
separated frequencies and comparing the results in view of the l/f 2 dependence 
of ionospheric group delay. When one undertakes such measurements the prac- 
tical problems of calibration, data extrapolation, time tagging and data com- 
pression become a primary concern. As usual, relative measurements are 
rather straightforward and absolute measurements become increasingly diffi- 
cult as higher accuracy is sought. 

An example of two frequency ranging is shown in Figure 14 where measure- 
ments were made at C-Band (4-6 GHz Applications Technology Satellite Ranging 
System) and at VHF (140 MHz). The tracking was between Rosman, N. C. and 
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the geostationary satellite ATS-3. The ionospheric diurnal effect is clearly 
traced out and the vertical total integrated content is calculated from the range 
difference data using equations 12 and 13. 


3.4 FARADAY ROTATION MEASUREMENTS 

Another widely used technique for observing integrated electron content is 
by means of Faraday rotation measurements. Once the ambiguity has been re- 
solved, a continuous record of total vertical integrated content is available. 
Again this value of total electron content accuracy is subject to certain much 
discussed qualifying approximations (Ref. 36, 37, 38 and 39). However the basic 
scheme is as follows. 

The Faraday rotation of the plane of polarization of a linearly polarized 
wave is given by the well known expression 


n = 


0. 0297 
f 2 


f 

•'n 


N(h)Mdh Radians 


(19) 


where M = H cos 6 sec i 

f = frequency of the wave 

N(h) = electron density at a height h in the ionosphere . 

H = the geomagnetic field intensity in MKS units 

8 = the angle between the ray and the geomagnetic field at the height h 
i = the angle between the ray and the vertical at the height h 
h s = vertical height of the satellite 

It is possible to approximate equation (19) by 


f 2 



N(h) dh 


Radians 


( 20 ) 


where M is the value of M at some ionospheric height chosen so as to make 
equations (19) and (20) numerically most equivalent. In many calculations the 
value of M is chosen to be that of M at a height near 400 km along the ray path. 
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A typical record of Q , the polarization rotation, is shown in Figure 15. 

This particular data was observed at GSFC on 25 February 1970 and during the 
solar eclipse of 7 March 1970 (Ref. 40). 

In this GSFC test Faraday rotation of the ATS-3 VHF signals was recorded 
by means of a two channel phase lock receiver, each channel having as an input 
one of the two circular components derived from the incoming linear polarization. 
The linearly polarized signal is resolved into right and left hand circular com- 
ponents and then multiplied in an electronic phase detector whose full scale output 
is 2n radians where each v of rotation in spatial angle results in one cycle of 
electrical phase output. The Smyth Research Associates (SRA-524) receiver 
output, as shown in Figure 16, was recorded directly by an Esterline Angus chart 
recorder running at a rate of 1.5 inches per hour. 

As shown in the records presented in Figure 15, during a normal day (such 
as February 25, 1970) ionization increases steadily during the morning hours up 
to the afternoon peak and then decays steadily as solar illumination decreases. 
Again with reference to Figure 15, a pronounced phase reversal occurred during 
the eclipse as the lunar shadow intercepted the line-of-sight between the GSFC 
site and ATS-3 and ionization abruptly decreased. At approximately 1400 hours 
local time the integrated electron content again increased up to the usual diurnal 
reversal. 

Figure 17 indicates the observed variation in total electron content which 
can be expected over any given month. This type of data clearly establishes the 
need for solar flux updating if accurate ionospheric modeling is desired. As 
indicated in Figure 4 this capability has been provided in the NASA-GSFC model. 
The dotted line in Figure 4 for measured Faraday rotation and/or f 0 F 2 input 
indicates an alternate option for improving worldwide predicted ionospheric 
effects. 

Figure 18 indicates the generally good agreement noted between modeled 
and measured integrated content extremes over a given month (Ref. 10). The 
measured data was obtained from Stanford as recorded in 1968 during the peak 
of the solar cycle. 

It should be pointed out that orbit computation typically involves the use of 
several days of tracking data with each day including coverage from several 
stations. Hence ionospheric corrections must be good "on the average". That 
is, the residual variations implied by, for example, Figure 18 are even further 
smoothed by the trajectory fitting process. 
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ELECTRON CONTENT (10 15 electrons/meter 1 ) 




NASA GSFC MODEL 

STANFORD FARADAY DATA 






LOCAL TIME (hours) 


3.5 TRAJECTORY ANALYSIS 


The range data is generally the most affected data type during any trans- 
ionospheric propagation. The bias is directly proportional to intervening electron 
content and independent of satellite motion. Therefore one method of viewing 
trajectory improvement is to calculate the RMS range residuals over a given 
data arc, comparing trajectories computed with and without ionospheric correc- 
tions. If VHF tracking is considered then the tropospheric range bias will be 
typically orders of magnitude less than the ionospheric bias. That is, the fre- 
quency independent tropospheric radio frequency bias at zenith is on the order 
of 2 to 3 meters and is also to a large extent modeled out (Ref. 42). Table 2 
shows the significantly improved range data residuals during two typical orbit 
computations for the AIMP-E and IMP-G satellites. Individual values of ion- 
ospheric range bias during specific tracking station coverage times will 
vary from 50 meters to several kilometers depending on time of day, geographical 
location, aspect angle and so on. 

Another means of detecting trajectory accuracy improvement is by means 
of trajectory overlap computations. Table 3 shows two seven day IMP-G satellite 
data arcs overlapped and compared again with and without the NASA-GSFC 


35 


a i i i a a a a u i u it ii ii a i I I i i 



Table 2 

Goddard Trajectory Determination System 
Test of NASA-GSFC Ionosphere Model 



VHF GRARR 
Tracking Date 


RMS Range Residuals 
(Meters) 

Spacecraft 

Length of Arc 

No Ionosphere 
Model 

Ionosphere 

Model 

AIMP-E 
(Lunar Orbit) 

Feb 25, 1972 

3 Days 

451 

140 

IMP-G 

Oct 18, 1971 

6 Days 

410 

200 


Table 3 

Orbit Accuracy Analysis 


Position Difference of Overlapping Arcs 


Maximum 

(km) 

Minimum 

(km) 

RMS 

(km) 

No Atmospheric 
Correction 

H 

0.4 

2.0 

Atmospheric 

Correction 

Applied 

2.5 

0.3 

0.7 


Note: 

1. Tracking of IMP-G at nominal 140 MHz 

2. Data Span 9 Nov. to 19 Nov. 1970 

3. Two seven day overlapping arcs compared 

4. Range and range rate tracking from 5 stations 
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Ionospheric Model. In this case the RMS of the absolute position difference 
relative to an Earth Centered Coordinate system was computed for the two over- 
lapped arcs. This measure of accuracy indicates an RMS error of 2.0 km without 
the model reduced to 0.7 km with the ionospheric model. The maximum difference 
was reduced from 7.4 km to 2.5 km. This is a marked improvement considering 
the fact that other modeling errors are present such as those associated with 
gravitational uncertainties, solar pressure uncertainty and station location errors. 


4.0 CONCLUDING REMARKS 

The NASA-GSFC Ionospheric Model for satellite tracking data corrections 
is part of the ongoing radio science and applications work being performed by 
the Geodynamics Program Division - GSFC to improve orbit computation ac- 
curacy. The NASA-GSFC ionospheric model as implemented in the operational 
support of the Radio Astronomy Explorer-B (RAE-B) has proven to be highly 
successful in improving data quality and hence trajectory accuracy. 

It cannot be overstressed that this work is simply an extension and practical 
application of work performed by numerous researchers beginning perhaps with 
the first sun spot observers in 1775! Although this ongoing NASA effort is pri- 
marily directed toward satellite tracking data improvement it also has direct 
application to such important topics as: 

• Correlation of ionospheric radio effects with satellite atmospheric 
drag modeling 

• Very Long Baseline Interferometry radio star refraction correction 

• Conventional radio astronomy refraction correction 

• Measurement of interplanetary plasma based on group-phase 
phenomena 

• Worldwide continuous predictions for 3 to 30 MHz skywave radio 
communications 

• Solar geophysical research 
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